Calculation of dephasing times in closed quantum dots 
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Dephasing of one-particle states in closed quantum dots is analyzed within the framework of 
random matrix theory and Master equation. Combination of this analysis with recent experiments 
on the magnetoconductance allows for the first time to evaluate the dephasing times of closed 
quantum dots. These dephasing times turn out to depend on the mean level spacing and to be 
significantly enhanced as compared with the case of open dots. Moreover, the experimental data 
available are consistent with the prediction that the dephasing of one-particle states in finite closed 
systems disappears at low enough energies and temperatures. 
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Quantum coherence of electrons in closed quantum 
dots has attracted much interest in recent years jj], 0, ||, ^, 
U . Electron-electron interactions are believed to be one 
of the main sources of dephasing in disordered systems at 
low temperatues. Compared to low-dimensional metals 
and semiconductors [|| , substantial modifications of this 
dephasing mechanism are caused by the confinement of 
the quantum dot which leads to discrete energy levels. 
In particular, the dephasing rate was predic ted P] to dis- 
appear at low excitation energies, e < Ay/ gj In g where 
A is the mean level spacing and g is the dimensionlcss 
conductance of the dot. 

Whereas there are a number of ways to measure the 
dephasing times in open quantum dots |7| |), the situ- 
ation is much more complicated in closed dots. Only a 
few experiments have attempted to study dephasing in 
closed quantum dots. Most of these have focused on the 
relaxation of highly excited states, verifying the continu- 
ous to discrete spectrum transition at e oc gA ||. Some 
signatures of dephasing in thcrmalized states have been 
studied by Patel et al. ||, who analyzed the statisti- 
cal distribution of the conductance maxima G max (the 
height of the Coulomb blockade peaks). They found 
that the ratio of standard deviation to mean peak height 
cr(G max )/(G max ) is smaller than what random matrix 
theory (RMT) predicts |To| , and attributed this reduc- 
tion to dephasing effects. More recently, Folk et al. Q| 
suggested to use the dependence of the conductance upon 
applying a magnetic field B, 

a = «G max >B#o - (G max ) s=0 )/<G max )^ , (1) 

as a probe of dephasing times. This is the closed dot ana- 
log of the weak localization magnetoconductance which 
was analyzed earlier for open dots 0. RMT predicts 
a = 1/4 Jll], [l2| , while Folk et al. measured consider- 
ably lower values of a, down to a w for the largest 
quantum dot with the maximal ratio fc^T/A (T is the 
temperature, ks the Boltzman constant), interactions. 
Beenakker et al. B theoretically analyzed the situation 



in which the phase-breaking inelastic relaxation rate r; n 
|l3f far exceeds the mean tunneling rate (inverse dwell 
time in the dot) T. It turns out that, in this limit, a is 
reduced much stronger than what found experimentally. 
Thus, they concluded that in the experiment [Q < T. 
However, as noted in Refs. [|| and || , the lack of a quan- 
titative theory of the crossover regime T ln ~ T, prevents 
a, full analysis of the experimental results. 

In this Letter, we study theoretically the effect of arbi- 
trary inelastic scattering on the conductance of a closed 
quantum dot. We develop an analytical approach that 
allows to evaluate a [Eq. (|l|)] and compare the results 
to the numerical solution. The approximate results are 
found to reasonably describe the behavior in the experi- 
mentally relevant temperature regime. Our calculations 
allow for the first time to extract dephasing times of 
low lying (thermaly excited) states in closed quantum 
dots. We observe a clear enhancement of the dephasing 
times relative to earlier results for open quantum dots . 
Moreover, contrary to the analysis of open quantum dots 
which showed a dependence on temperature alone, we 
find a dependence on both T and A. From our analysis it 
follows that the measurements of Folk et al. B arc not 
inconsistent with vanishing dephasing rate for low exci- 
tation energies Q . A more detailed presentation of the 
calculation will be given in JlJ] . 

In the experimentally relevant regime KTi n , HT < 
ksT, A, each state of the quantum dot is determined 
by a tuple {rii} of occupation numbers for the single 
particle eigenstates with energies Ei and spins Si. The 
electrons can tunnel between the dot and the two leads. 
The left (L) and right (R) leads differ due to the ap- 
plied voltage V. The elctrons in the leads are ther- 
malized and distributed according to the Fermi function 
Ifd{E) = [l + e^^]- 1 . The probability PaKK}) to 
find a given set of occupation numbers {rii} with a total 
of M electrons (restricted to M 6 {N, N + 1} due to the 
Coulomb blockade) obeys the following Master equation 



2 



dP N+1 ({rn}) 
dt 



=Y^Sn i oT$[(l-ti)P N+1 ({n i }^)-tiP N ({n i })} +£^ W [itiMW}^) - rg> w ({n,})] (2) 



J A 



r 



Here, {ni} + j ({nj}_j) are the tuples obtained from {rii} 
by adding (removing) one electron in the one-particle 
eigenstate j, and ff = fFD(Ej + {5 X l - l/2)eV - At), 
where A € {L,R} and fj, is the effective chemical poten- 
tial, including the charging energy. 

The first terms in (|J) describe the tunneling of elec- 
trons between the dot and the leads. The additional 
terms F^ in (||) describe inelastic scattering processes 
between the dot's one-particle eigenstates j and k. We 
assume that these are caused by thermal bosonic fluctua- 
tions at temperature T and neglect any back-coupling of 
the scattering to this bose bath. Under the assumptions 
that the coupling strength is independent of the specific 
levels involved and there is no spin-scattering, one arrives 
at (uijk = Ej — Ef.) 
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(3) 



m exp[u jk /(k B T)] 

where D(E) is the bosonic density of states. As will 
be shown below, the suppression of a is quite robust to 
the specific model of interaction, and depends mainly 
on the total inelastic scattering rate Ti n . We consider 
Tin as a phenomenological parameter, to be determined 
experimentally. Since the experimental quantum dots 
j| contain a large number of electrons, N 3> 1, they 
can be described by RMT jf6). In particular, the tun- 
neling rates are Porter Thomas distributed Pp(T) cx 



exp[-/?r/(2r)], with 0=1,2 for the Gaussian or- 



thogonal (no magnetic field) and unitary ensemble (with 
a magnetic field). The difference between P\(r) and 
P2(r) leads to the afore mentioned of a = 1/4 in the 
absence of inelastic scattering 
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The inelastic scattering model (|3|) is exponentially cut 
off to states outside an energy window of 0(k B T) and, 
thus, Tin vanishes at low temperatures. At fc^T ^> A 
on the other hand, there are many states M oc T/A 
connected by the inelastic scattering. Therefore, for T — > 
co, the total inelastic scattering rate r; n /r — > oo and the 
result of || is approached. In the leading order in A/ k B T 
the solution of the Master equation (0) is the equilibrium 
distribution and 
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i.e., a = (here, T is the mean tunneling rate to lead A, 
r = V + r ). Two corrections arise in the next order in 



A/k B T: (i) the effect of a finite total inelastic scattering 
rate and (ii) T is replaced by an average M levels around 

the Fermi energy, i.e., f A = £^(1 - f 3 )^ 

H |l5) (Here and in the following fj is the Fermi func- 
tion at both leads for V = 0). This introduces correc- 
tions 0(1/ M) due to correlations between the nomina- 
tor and denominator of Eq. (Q). We calculate the for- 
mer and take into account the latter by solving the Mas- 
ter equation using perturbation theory in r/r; n where 
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where P eq (N) is the equilibrium probability to have N 
electrons in the dot and r* n = r^ n /(l — fj). We neglect 
the weak j-dependence of r* n . The total inelastic scat- 
tering rate T m is defined as value of Tj n at Fermi energy, 



i.e. 



r? for Ei 



0. 



In the following, we employ an approximation inspired 
by this high-T expansion. To first order in voltage, one 
can write Pjv({n,}) = P^ q ({nJ)(l + ^ £\ <W*(i)), 
where Pt/({ni}) is the equilibrium value of Pj\f({ n i})- 
With a large number of final states to scatter to, one can 
replace the sum over many individual terms P({nj}+fc_j) 
in the Master equation (||) by an averaged quantity 

Y J ^iPN + x{{n i } +k -j)-V^P N+l {{n i }) 



^i(W)£l£[¥(fc) -*(?)] 



eV 



(7) 



Here, ^ should, in principle, be a weighted average over 
levels within a range of 0(T) around a particular level 
j considered. However, only levels around the Fermi en- 
ergy are of interest for the conductance since the contri- 
bution of every level j to the conductance is multiplied 
by fj(l — fj). For this reason, we approximately employ 
a constant * = £\ f 3 (l- f 3 )^U)/(E J /iCWi)) inEqs. 



3 




FIG. 1: Comparison of the numerical solution of the full Mas- 
ter equation with the high temperature approximation. The 
lattter is seen to work well for ksT > A. 

and (^) This leads to a self-consistent solution of the 
Master equation with the result [jl4| 

a (■?+ (ffiffir >)}■> 

where r* ot = (rf + Tf + r* tt ) _1 . One would obtain the 
same form (||) but with ((...)) = J2jLi ■■■ considering M 
degenerate levels filled with N e {0, 1} electrons. 

The result (||) can be interpreted in the following way: 
the first term represents processes in which the electron 
was not scattered at all. These happen with probabil- 
ity (Tf + rf )r' ot and the resulting conductance peak 

pi pi? 

heights are proportional to p-fe^ ; yielding rfT f"T* ot al- 
together. The second term represents contributions from 
electrons that were inelastically scattered after tunneling 
from one lead, and their contribution to conductance is 

Equation (Q) is the main result of this paper. It is 
based on an approximation (^) which can be justified in 
the high temperature limit. The particular advantage of 
this approach is that it gives not only the correct leading 
high temperature behavior [Eq. (||)] but also reproduces 
correctly the limits r; n = and Ti n = oo for all T includ- 
ing a = 1/4 at T = 0. Below we demonstrate that this 
approach works pretty well in the intermediate regime 
k B T ~ A. 

In order to calculate G and a one has to average 
Eq. (||) w.r.t. the different ensembles. One could do 
so numerically but it is possible to get analytical re- 
sults via expanding Eq. (||) in powers of A/ksT Q. 
The first three terms in the A/ksT expansion already 
give good accuracy in the relevant regime fcsT > A 
and are employed in the following. As we are inter- 



FIG. 2: A contour plot of a as a function of T/A and Tin/r, 
based on the high temperature approximation. The values 
the bold contours are specified. Given T, A and a from future 
experiments, one can extract Fi n /T from this Figure. 

ested in this regime we assumed a picket fence distribu- 
tion with spacing A between consecutive spin-degenerate 

levels (E 2j = E 2j ^ = jA; I* = I* _ l5 r* = f/2). 

We tested the range of validity of this high- 
temperature approximation against the numerical solu- 
tion Jl8| of the Master equations (|J). The latter is ob- 
tained by solving the Master equation (||) by sparse ma- 
trix inversion [ fl8| . Figure ^ compares values of a, as 
calculated using the first three terms in the A/T expan- 
sion, with the numerical values. The agreement is very 
good for sufficiently high temperatures, and reasonable 
even for low T. In the whole temperature regime, the 
deviations are within current experimental accuracy. 

It, thus, appears that our analytical approach provides 
a reliable way to determine Ti n from the experimental 
measurements of a, in the whole temperature regime. 
For future experiments we provide Fig. |^, which presents 
a as a contour-plot in the space spanned by fcsT/A and 

r in /r. 

A direct experimental test would be provided by mea- 
suring values of a in a given dot at fixed T, as a function 
of r (which can be achieved by changing the contact set- 
ting). The theoretical dependence of a on T involves a 
single fitting parameter, i.e., the unknown total scatter- 
ing rate Ti n which is assumed to be unaffected by the 
contact setting. A first step in this direction was done in 
Q , and in the inset of Fig ^ we compare the prediction of 
our high temperature approximation with the measure- 
ments of a for three different values of T. An excellent 
agreement is obtained, though more data points are re- 
quired for reliable conclusions. 

We now use the above theory to extract dephasing 
times from the data points (mean values and error bars) 
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FIG. 3: Dephasing times, t^, as extracted from the data 
points in Ref. [^] for four different dots: A = 28/xeV (cir- 
cles, long-dashes error bars), A = 10/xeV (squares, solid er- 
ror bars), A = 2.4/ieV (up-triangles, dashed error bar), and 
A = 0.9/xeV (dot-dashed error bar); dotted line: fit so open 
dot experiments as calculated in |7J . Error bars which extend 
up (down) beyond the graph should be understood as going 
up to infinity (down to zero); if no corresponding point is 
visible the experimental mean value itself gives = oo (or 
T,f, — 0). In the inset, we fit experimental measurements for 
different values of P Q with our theory. Fhe single fitting 
parameter is Fi n = 0.25^ieV, or T0 = 16ns. 



of Folk et al. Q. Figure || presents these estimates as 
symbols and error bars, respectively, and compares them 
with open dot values |Q. A clear enhancement of the 
dephasing times compared to open dots is observed. In 
addition, dephasing times strongly depend on A (as can 
be seen at T = 45 mK). This is in contrast to open dot 
results Jt|. An additional suppression of a for fcsT < A, 
resulting from level-spacing fluctuations [lTj was not 
included in our analysis, and therefore our results under- 
estimates the dephasing times for ksT < A. In addi- 
tion, the result for the A = 0.9/ieV quantum dot which 
is consistent with = should be interepreted care- 
fully since the result implies HTi n > A and the Master 
equation is not applicable anymore. Based on our anal- 
ysis, t4ifireceirtra 



on measurements of the weak-localization correction a. 
Analyzing a recent experiment by Folk et al. 0), we 
see a clear enhancement of the dephasing time compared 
with open dots values. There is no inconsistency with 
theoretical predictions for electron-electron interaction, 
in particular, a vanishing dephasing rate at a critical 
A — dependent temperature. We note, however, that 
the available experimental data is limited and has con- 
siderably statistical uncertainties. Future experiments 
are necessary and we offer Figure || to extract the tem- 
perature and level-spacing dependence of the inelastic 
scattering rate and to thoroughly test the prediction of 
a diverging dephasing time. 
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closod quantum dots is consistent with dephasing due to 
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electron-electron interaction alone, including the predic- 
tion of the critical vanishing of dephasing rate. However, 
given the large error bars of the current experimental 
data, one can not exclude an algebraic behavior or even 
a saturation of the dephasing rates for T — > 0. Neverthe- 
less, the behavior is clearly different from that of open 
quantum dots (7) and is A-dependent. 

In conclusion, we provide a theoretical approach to ex- 
tract the inelastic scattering rate in closed dots based 
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[17] In the temperature regime 0.1 A < T < A, differences 
in the GOE/GUE level spacing distributions lead to a 
significant reduction from the value a = 1/4 (down to 
a « 0.2) ^j. This effect alone can explain the suppres- 
sion of a in the experiment 0] for ksT < A. 
[18] For ksT < 1.6A, we took into account all configura- 
tions {ni} which involve levels in [-4fcsT..4fcsT] with 
P e<l {ni}/mm {ni} P cq {n,} > exp(-5A/fc s T); for k B T > 
1.6A the interval was [— 3fcsT .. 3fcsT] and the exponen- 
tial cut-off -4.5A/fc s T, see Q for details. 



